New Bardeen- Cooper- Schrieffer- type theory at finite temperature 

with particle-number conservation 



H. Nakade 

Department of Physics, Faculty of Science, 
Chiba University, Inage, Chiba 263-8522, Japan 

K. Tanab^ 

Department of Physics, Faculty of Science, 
Saitama University, Sakura, Saitama 338-8570, Japan 
(Dated: February 1, 2008) 

Abstract 

We formulate a new Bardeen-Cooper-Schrieffer (BCS)-type theory at finite temperature, by 
deriving a set of variational equations of the free energy after the particle-number projection. 
With its broad applicability, this theory can be a useful tool for investigating the pairing phase 
transition in finite systems with the particle-number conservation. This theory provides effects of 
the symmetry-restoring fluctuation (SRF) for the pairing phenomena in finite fermionic systems, 
distinctively from those of additional quantum fluctuations. It is shown by numerical calculations 
that the phase transition is compatible with the conservation in this theory, and that the SRF shifts 
up the critical temperature (T cr ). This shift of T cr occurs due to reduction of degrees-of-freedom 
in canonical ensembles, and decreases only slowly as the particle-number increases (or as the level 
spacing narrows), in contrast to the conventional BCS theory. 

PACS numbers: 05.70.Fh, 21.60.-n, 74.20.Fg, 05.30.Fk 
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Pairing phenomena have been observed in various fermionic systems. Condensate of 
fermion pairs is realized at low temperature (T), and is usually described by the Bardeen- 
Cooper-Schrieffer (BCS) theory 0] or the Hartree-Fock-Bogolyubov (HFB) theory As T 
increases, the pairing phase transition occurs; the condensate is dissolved at a critical tem- 
perature T cr . In the pairing theories such as BCS or HFB, the particle-number (n) conser- 
vation is violated in the low T (i.e. superfluid or superconducting) phase, as a result of the 
spontaneous symmetry breaking. While this picture is reasonable if n is practically infinite, 
certain interest has been attracted by finite systems such as atomic nuclei Q, in which the 
condensate of fermion pairs forms at low T but the n conservation is actually preserved. We 
have examples also in mesoscopic systems; e.g. granular films [3|, superconducting islands 0] 
and ultrasmall metallic grains It has been argued that there are no sharp transitions 

in finite systems, whose signatures are washed out by quantum fluctuations jj • Efforts have 
been made to find fingerprints of the superfluid-to-normal phase transition in atomic nuclei. 
Recent experiments have revealed an S-shape in the graph of the heat capacity C as a func- 
tion of T, which was extracted from high-precision level-density measurements jjj. While 
this S'-shape was suggested to be the fingerprint, we have shown that the n conservation 
produces a similar S'-shape even without the transition JjJ. Under this situation the follow- 
ing fundamental questions are raised for phase transitions and their relation to finiteness: 
(i) for increasing T what roles the conservation law plays in finite systems, and (ii) how a 
'phase transition' develops as n increases. To answer these questions, it is desired to apply 
the n projection in the pairing theory, particularly in the variation-after-projection (VAP) 
scheme. It was shown that exact results in canonical ensembles (CE) are well approximated 
by the n-projected BCS approach in the VAP scheme, for the degenerate model with a con- 
stant pairing jlp| . In this Communication, we formulate a new BCS-type theory at finite 
temperature with the n projection, which is well founded on the variational principle and 
has much wider applicability. 

The BCS theory is closely linked to the Bogolyubov transformation c\. = UkCi\ — VkCij.. 
Here c k stands for the creation operator of the original fermion on the single-particle (s.p.) 
state k, a, the creation (annihilation) operator of quasiparticle (q.p.) on k, and k 

represents the time-reversal to k. The unitarity derives Uk = \/l — \vk\ 2 , and Vk is usually 
taken to be real. At finite T, the BCS theory is obtained for grand-canonical ensembles 
(GCE), by assuming the trial statistical operator 11] 

WG = Tr(e- g «/ T ) ' H ° = 2^ £ k a k<Xk, (1) 

where we set the Boltzmann constant k-Q = 1 and the trace is taken over the GCE. This is an 
approximation on the excitation spectra of the system by those of the non-interacting q.p.'s. 
The parameters Vk and or equivalently fk = l/(e £fe//T + 1), are determined by variation 
of the free energy. This theory is specifically called GCE-BCS theory in this paper. 

To include a certain part of n conservation effects, the number-parity (7T n ) projection [l^ 



2 



has conveniently been applied. However, formulation of the full n projection, by which we 
can work in CE, has been limited for the variation-before-projection (VBP) scheme (lH 
Q], or for the degenerate model with a constant pairing except at zero T. Whereas 
application of the VAP scheme is desired to investigate effects of a conservation law as stated 
earlier, it has been difficult because of a problem generic to VAP treatment of spontaneously 
broken symmetries. We exemplify the problem for the n projection case. The n projection 
operator is written as 

P n = i e** N ^>dip , (2) 
^ Jo 

where N = J2k-^k, N k = c\ct- With the approximate trial statistical operator 

p p-Ho/Tp 

WC ~ Tr(e-^/T Pn) ' ( 3 ) 
the free energy in the BCS-like picture can be defined by 

F c = E c -TSc; E c = Tr( Wc #) , S c = -Tr(w c In w c ) . (4) 

However, Sq is not tractable in general cases, because [P n ,H ] ^ and lnu> c leads to an 
infinite series of P n e~ H °l T P n . 

To derive a VAP equation in CE, it is practical to introduce an additional approximation 
on the entropy 

F C = E C -TS C ; S C = ^XvpS 1 +^- H °> Tp «) ■ («) 

The Peierls inequality 0| Fq > Fq > _F exact holds, which justifies variation of Fq. Although 
Sq can be negative at T m 0, violating the third law of thermodynamics, the variation may 
suppress influence of this problem. 
We here use the following expression, 

M, = Tar; d* = ^- n) (IIcn^. ( 6 ) 

J k>0 

with C = fk(l - h) + (1 - fk)h + fkMT + (1 - /*)(! - er = «2e^ + % 2 e-^, and 

Q denotes half the number of the s.p. states. Although Q can be infinity in principle, we 
should take it to be finite in practical applications, by introducing a proper cut-off. The 
expression k > for FJ indicates that the product is taken once for the (kk) pairs. We then 
have, for a given n and T, 

™ c = + r£[^%, fc , (7) 

k k>0 k 

Me = £[/« + /i£^ (8) 

k>0 k>0 
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where / fc c = Tr(w c aja fc ), = Ti[w v H) with w* = e~^ N r e~ H °/ T ^/Tr(e"^ 7V r e~ H °/ T ), pi 



Tr (w*N k ), kI = Tr(w^c k c k ), k\ = Tr^cjcjh, and 



K 
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-A* = ^ 
k SkI ' 



-A*=*^ 



(9) 



Owing to the extended Wick's theorem [3], E v can be expressed in terms of p£, kI and k\°. 
The variation of F c with respect to f k leads to the following coupled equations of e k , 

E £ *'§¥ = [W-(^ + ^)(A-(A-«*)^)+^Ajf e -^ + ASf e ^}/^] 

-[(^-^ + ^c)^],, (10) 

where D* = Kp\ + h\p\ - A£*jf - Aj>*£. Since / fe c = / fc + / fc (l - /*) [^]^, for the lhs 
of Eq. (JTDJ) we have 
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1 — COS (p 



<W { 1 + (1 - 2f k ) ~ M-wfrO- ~ /fc) [ 



The variation of Fc with respect to v k yields 



2u k v k h k - (u\ - vl)A k = ; 
Afc = ^(A^ + A^)}/Crl 

-itp 



(12) 
(13) 



hk = [^ht + ht)^ + (D*-E* + E c + T)^] v 



/^Ek'fk'jl ~ /fc') 



#fcfc' + <5fcfc' 
1 — /fc — /fc 



— (i — <W — <W) 



5/fc' L 9/fc/ V 
% sin 



x- 



Sfc 



(14) 



Equation (fT^j) can be solved in analogy to the usual BCS equation, 
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hi + A? 



(15) 



In this regard A k plays a similar role to the gap parameter in the usual BCS theory. 

We shall call the above formalism canonical- ensemble BCS (CE-BCS) theory. The CE- 
BCS theory recovers the n-projected BCS theory for the ground state in the T — > limit, and 
coincides the GCE-BCS theory when g?$ is replaced by S(ip) dip. From the GCE-BCS theory 
viewpoint, the (p variable in Eq. (j2J) corresponds to the Nambu-Goldstone mode which the 
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broken symmetry gives rise to. The (p integration takes account of the quantum fluctuation 
of the Nambu-Goldstone field. Since the n conservation is restored by this fluctuation, we 
call this fluctuation symmetry-restoring fluctuation (SRF). Via the n projection the SRF is 
separated from additional quantum fluctuations (AQF). The SRF may change configuration 
{e.g. Vk) that minimizes the free energy. This effect is taken into account in the CE-BCS 
framework, not in the VBP scheme. It is also noted that the CE-BCS theory is broadly 
applicable with no limitation on Hamiltonian, and is relatively easy to increase n or to 
enlarge the model space, compared with other extensions of the BCS theory that restore the 
n conservation. 

We present a numerical application of the CE-BCS theory. The following Hamiltonian is 
adopted for the sake of simplicity, 



In any of the examples below, g is adjusted so as for the gap parameter of the GCE-BCS 
approximation (which is denoted by Aq) to be unity at zero T. In other words, all quantities 
having the energy dimension are represented in unit of Aq(T = 0). The model space is cut 
off by \tk\ < A, and we set A = 10. The parameter is chosen to be tk = —A + (k — l)d 
where d = 2A/ (fl — 1) (k = 1,2, ■ ■ ■ with the time-reversal symmetry tj. = tk- Keeping 
n = Q (i.e. half-filled), the particle number n is varied. In this paper we restrict ourselves to 
n = even cases. The parameter /i in Eq. (JTfij) merely shifts the zero-point of energy in CE, 
while it corresponds to the particle number condition in GCE. Indeed, though not trivial 
in Eqs. (jl0H4|) . the CE-BCS results hardly depend on /i. For fast and safe convergence in 
numerical calculations, minimization of Fq is implemented by combining Eqs. (|10I15J) with 
the steepest decent method. 

The approximation on the entropy by Sc has been tested for the n = 26 case. Since 
it is difficult to compute Sc, we compare Sc with the exact canonical entropy S exact = 
— Tr(uTnu>) where w = P n e~ H ' T P„ /Tr(e~ H / T P n ) , which can be obtained b y t he quantum 
Monte Carlo (QMC) calculation [17[, using the method described in Ref. [18]. Although 
there is a slight discrepancy at low T as mentioned earlier and a weak kink around T = T C T 
(Tq will be defined later), Sc is found to be in moderate agreement with S exact at any T. 

The thermal expectation value of an operator O in the CE-BCS (GCE-BCS) is denoted by 
(O)c ((O)g). Replacing P n by the 7r n projection operator [14], we can calculate 7r n -projected 
counterparts to the n-projected quantities. The 7r n -projected expectation value of O will be 
considered for comparison, which is expressed as {0) n , as well as the expectation value in the 
VBP (0)c'(= {OP n )o). The gap parameter A G (= g(B) G ) is regarded as an order parameter 
for the pairing transition in the GCE-BCS theory. However, obviously (B)c = (B)c — at 
any T. Instead we consider two alternative definitions; one is A^ in Eq. (113)1 . and the other 
is A av = g^{BlB) - }2 k >o( N k) ( N k) B- For the latter, we define A c v , A$ and so forth, in 
accordance with the expression for ( ) . These pairing parameters have lost direct connection 




(16) 
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to the energy gap. In Fig. E^a), Aq(T) and Aq(T) are depicted for various n values. For 
comparison, Aq,(T) and A^(T) are presented for the n = 26 case. Because A& depends 
on k, A min (T) = min/c (A^(T)) and A max (T) = max^ (Afc(T)) are shown for n = 26. It is 
confirmed that A^'s and Aq v do not differ much from one another. 

The gap parameter Aq vanishes at the critical temperature Tq(^ 0.6) in the GCE-BCS 
theory. As n increases Aq converges rapidly; Aq's for n = 26 and 56 cannot be distin- 
guished in Fig. HJa). It is remarked that, in the CE-BCS results, transition-like behavior 
remains. A^ v vanishes at a certain temperature Tq, which is appreciably higher than Tq. 
Tq approaches Tq as n grows, but only gradually. Although shift of the critical temperature 
is already found in A^ v (giving T™), it is far less significant than in A^ v . While the SRF 
enhances the pairing parameter in the CE-BCS results, the VBP scheme gives smaller Aq, 
than Aq. Effects of the SRF are not carried correctly in the VBP, since the configuration 
change is discarded. In Fig. W[b), we compare the (B'B) values in the BCS-type approx- 
imations with the exact ones for n = 26, which are evaluated by the QMC calculation. 
Taking account of the SRF effects, (B^B)qs are close to the exact values at T < Tq, unlike 
(B^B)q, (B^B)c' and (B^B) n . Thus the CE-BCS theory gives significant improvement over 
the GCE-BCS and the 7r n -projected theories, while keeping simplicity of the BCS picture to 
a considerable extent. The deviation at T > Tq may be attributed mainly to influence of 
the AQF. 

The above behavior of the pairing parameters is reflected in the q.p. energies. Although 
e/c's are parameters independent of v^s, they also correlate well to the pairing parameters, 
having a certain energy gap at T < Tq, as shown in Fig. life). Since is expected to give 
an approximate energy of the adjacent odd-n system, this suggests even-odd difference in 
energy at T < Tq. 

We next view the heat capacity C(T) = dE/dT (C G for GCE-BCS, C c for CE-BCS, 
and so forth), whose singular structure, if any, is linked to a phase transition in general. 
Figure |2fa) shows a specific heat C(T)/n, where C is computed by numerical differentiation 
of E = (H). As Cq at Tq, Cq has discontinuity at Tq. It is thus fair to say that the pairing 
transition remains in the CE-BCS approximation, but at Tq that is substantially higher 
than Tq. The SRF shifts up T cr , not erasing the transition, in the CE-BCS theory. The 
S-shape behavior of C(T) at low T may be compared to those observed & Comparison 
with the exact QMC result draws consistent consequence with the {B^B) case in Fig. IHb). 

In Fig. Efb), we present expectation values of the q.p. number Af = J2k a l a ki as a 
function of T. It is found that the shift from Tq to Tq is closely connected to the difference 
between (M)g an d (N)c- The n projection reduces the number of states. At low T, lack of 
one q.p. states gives rise to reduction of (A/ - )^ and (A/")c, in comparison with (A/")g M- At 
higher T, the n projection eliminates some of the higher q.p. degrees-of-freedom. Therefore, 
(A/}c becomes smaller than (J\f) n . Moreover, as T becomes higher, the entropy increases 
more slowly in CE than in GCE and in the 7r n -projected space. This further delays rise of 
(A/)c through the configuration change. This effect is the stronger for the fewer n, and leads 
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FIG. 1: (a) n- and T-dependence of pairing parameters. Colors distinguish n; green for n = 10, 
blue for n = 26 and red for n = 56. For each color the thick solid line represents A^, while the 
dashed line Ag- For the n = 26 case, A m i n and A max (dotted lines), Afy (dot-dashed line) and 
A| v (thin solid line) are also presented, (b) (B>B) for n = 26. Diamonds are the exact quantum 
Monte Carlo results. We use the same conventions for lines as in (a), (c) Ek (k — n/2 = 1,2,3,4) 
for n = 26. Conventions are the same as in (a). 
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FIG. 2: n- and T-dependence of (a) specific heat C/n, and (b) (N)/n. See Fig. ^for conventions, 
to the shift of T cr . 

Figure El depicts the n-dependence of T cr . Tq rapidly goes to the bulk (i.e. n — > oo) limit 
T™. On the contrary Tq approaches slowly. We find, by fitting, Tq—T^ oc n~ a75 oc d ' 75 , 
except for quite small n. For a fixed s.p. level spacing d, Tq is insensitive to the cut-off A, 
although it slightly goes up as A increases. 

The 'phase transition' picture originates from the approximation of H by Hq in the trial 
statistical operator. While the transition will be washed out due to the AQF, the present 
study yields the approximate picture in which the phase transition is compatible with the 
conservation law. The finite system approaches the bulk limit as both the SRF and the 
AQF are reduced, and the reduction of the SRF gives rise to the decrease of Tq. It should 
be commented, however, that in the degenerate model the exact treatment of Sc was 
shown to wash out the signatures of the transition, while the CE-BCS theory keeps the 
transition behavior due to the approximation of Sc by Sc- Even slight difference in the 
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FIG. 3: n-dependence of T cr ; Tq (pluses), (triangles) and (crosses). The solid line is a fit 
to Tg r ; T cr = T£ + 6.8 n~ - 75 . 



entropy could influence occurrence of the transition. Whereas the present study gives a new 
insight, further investigation is necessary for full understanding of the superfluid-to- normal 
transition under the n conservation. 

In summary, we have formulated a BCS-type theory in canonical ensembles (CE-BCS 
theory) at finite temperature, by introducing an approximation on the entropy with retaining 
the variational principle. The CE-BCS theory keeps some of the structure in the conventional 
grand-canonical ensemble BCS (GCE-BCS) theory. The equations in the theory are well 
connected to the GCE-BCS theory, so that, under the particle-number conservation, the 
GCE-BCS theory could be regarded as an approximation of the CE-BCS theory. Having 
broad applicability to a moderately good accuracy, the CE-BCS theory provides us with 
a useful tool to investigate effects of the symmetry- restoring fluctuation (SRF), separately 
from those of the additional quantum fluctuations, although a part of the SRF might be 
missed due to the approximation on the entropy. Numerical application of this theory gives 
a new picture for the pairing transition in finite systems, in which the pairing transition 
is reconciled with the particle-number conservation. In the CE-BCS framework the SRF 
has been found to shift up the critical temperature T cr . The shift of T cr occurs due to 
the reduction of excitation degrees-of-freedom in canonical ensembles. T cr in the CE-BCS 
theory approaches its bulk limit much more slowly than that in the GCE-BCS theory, for 
narrowing single-particle level spacing. Several more aspects of the CE-BCS theory that give 
us significant insights into the pairing phenomena will be discussed in a future publication. 
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